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ABSTRACT 

We develop a simple yet comprehensive method to distinguish the underlying drivers 
of galaxy quenching, using the clustering and galaxy-galaxy leasing of red and blue 
galaxies in SDSS. Building on the iHOD framework developed by Zu & Mandelbaum 
(2015), we consider two quenching scenarios: 1) a “/loZo”quenching model in which halo 
mass is the sole driver for turning off star formation in both centrals and satellites; 
and 2) a quenching model in which the quenched fraction of galaxies depends 

on their stellar mass while the satellite quenching has an extra dependence on halo 
mass. The two best-fit models describe the red galaxy clustering and lensing equally 
well, but halo quenching provides signihcantly better fits to the blue galaxies above 
IO^'/j^^Mq. The halo quenching model also correctly predicts the average halo mass of 
the red and blue centrals, showing excellent agreement with the direct weak lensing 
measurements of locally brightest galaxies. Models in which quenching is not tied to 
halo mass, including an age-matching model in which galaxy colour depends on halo 
age at fixed M*, fail to reproduce the observed halo mass for massive blue centrals. We 
find similar critical halo masses responsible for the quenching of centrals and satellites 
(~ 1.5 X IO'^/z^'Mq), hinting at a uniform quenching mechanism for both, e.g., the 
virial shock-heating of infalling gas. The success of the iHOD halo quenching model 
provides strong evidence that the physical mechanism that quenches star formation 
in galaxies is tied principally to the masses of their dark matter halos rather than the 
properties of their stellar components. 

Key words: cosmology: observations — cosmology: large-scale structure of Uni¬ 
verse — galaxies: luminosity function, mass function — gravitational lensing: weak — 
methods: statistical 


1 INTRODUCTION 

The quenching of galaxies, namely, the relatively abrupt 
shutdown of star formation activities, gives rise to two 
distinctive populations of quiescent and active galaxies, 
most notably manifested in the strong bimodality of galaxy 
colours (Strateva et al. 2001; Baldry et al. 2006). The under¬ 
lying driver of quenching, whether it be stellar mass, halo 
mass, or environment, should produce an equally distinct 
split in the spatial clustering and weak gravitational lensing 
between the red and blue galaxies. Recently, Zu & Mandel¬ 
baum (2015, hereafter Paper I) developed a powerful sta¬ 
tistical framework, called the iHOD model, to interpret the 
spatial clustering (i.e., the projected galaxy autocorrelation 
function Wp) and the galaxy-galaxy (g-g) lensing (i.e., the 
projected surface density contrast AX) of the overall galaxy 
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population in the Sloan Digital Sky Survey (SDSS; York 
et al. 2000), while establishing a robust mapping between 
the observed distribution of stellar mass to that of the un¬ 
derlying dark matter halos. In this paper, by introducing two 
empirically-motivated and physically-meaningful quenching 
models within iHOD, we hope to robustly identify the dom¬ 
inant driver of galaxy quenching, while providing a self- 
consistent framework to explain the bimodality in the spatial 
distribution of galaxies. 

Galaxies cease to form new stars and become quenched 
when there is no cold gas. Any physical process responsible 
for quenching has to operate in one of three following modes: 

1) it heats up the gas to high temperatures and stops hot 
gas from cooling efficiently (e.g., gravitational collapse and 
various baryonic feedback; see Benson 2010, for a review); 

2) it depletes the cold gas reservoir via secular stellar mass 
growth or sudden removal by external forces (e.g., tidal and 
ram pressure; Gunn & Gott 1972); and 3) it turns off gas 
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supply by slowly shutting down accretion (e.g., strangula¬ 
tion; Balogh & Morris 2000). However, due to the enormous 
complexity in the formation history of individual galaxies, 
multiple quenching modes may play a role in the history of 
quiescent galaxies. Therefore, it is more promising to focus 
on the underlying physical driver of the average quenching 
process, which is eventually tied to either the dark mat¬ 
ter mass of the host halos, the galaxy stellar mass, or the 
small/large-scale environment density that the galaxies re¬ 
side in, hence the so-called “halo”, “stellar mass”, and “envi¬ 
ronment” quenching mechanisms, respectively. 

Halo quenching has provided one of the most coher¬ 
ent quenching scenarios from the theoretical perspective. 
In halos above some critical mass (Mshock~10*^/i“*Mo), virial 
shocks heat gas inflows from the intergalactic medium, pre¬ 
venting the accreted gas from directly fueling star forma¬ 
tion (Binney 1977; Birnboim & Dekel 2003; Katz et al. 2003; 
Binney 2004; Keres et al. 2005, 2009). Additional heating 
from, e.g., the active galactic nuclei (AGNs) then maintains 
the gas coronae at high temperature (Croton et al. 2006). 
For halos with Mi,<Mshock, the incoming gas is never heated 
to the virial temperature due to rapid post-shock cooling, 
therefore penetrating the virial boundary into inner halos as 
cold flows. This picture, featuring a sharp switch from the 
efficient stellar mass buildup via filamentary cold flow into 
low mass halos, to the halt of star formation due to quasi- 
spherical hot-mode accretion in halos above natu¬ 

rally explains the colour bimodality, particularly the paucity 
of galaxies transitioning from blue, star-forming galaxies to 
the red sequence of quiescent galaxies (Cattaneo et al. 2006; 
Dekel & Birnboim 2006). To first order, halo quenching does 
not discriminate between centrals and satellites, as both are 
immersed in the same hot gas coronae that inhibits star 
formation. However, since the satellites generally lived in 
lower mass halos before their accretion and may have re¬ 
tained some cold gas after accretion, the dependence of satel¬ 
lite quenching on halo mass should have a softer transition 
across Mshock, unless the quenching by hot halos is instanta¬ 
neous. 

Observationally, by studying the dependence of the red 
galaxy fraction on stellar mass M, and galaxy environ¬ 
ment 55^^ (i.e., using distance to the 5th nearest neighbour) 
in both the Sloan Digital Sky Survey (SDSS) and zCOS- 
MOS, Peng et al. (2010, hereafter PIO) found that can 
be empirically described by the product of two independent 
trends with M, and ^jnn, suggesting that stellar mass and en¬ 
vironment quenching are at play. By using a group catalogue 
constructed from the SDSS spectroscopic sample, Peng et al. 
(2012) further argued that, while the stellar mass quench¬ 
ing is ubiquitous in both centrals and satellites, environment 
quenching mainly applies to the satellite galaxies. 

However, despite the empirically robust trends revealed 
in PIO, the interpretations for both the stellar mass and 
environment trends are obscured by the complex relation 
between the two observables and other physical quantities. 
In particular, since the observed M, of central galaxies is 
tightly correlated with halo mass Mi, (with a scatter ~0.22 
dex; see Paper I), a stellar mass trend of is almost in¬ 
distinguishable with an underlying trend with halo mass. 
By examining the inter-relation among M,, Mi,, and ^jnn, 
Woo et al. (2013) found that the quenched fraction is more 
strongly correlated with Mj, at fixed M, than with M, at Mi,, 


and the satellite quenching by (Isnn can be re-interpreted as 
halo quenching by taking into account the dependence of 
quenched fraction on the distances to the halo centres. The 
halo quenching interpretation of the stellar and environment 
quenching trends is further demonstrated by Gabor & Dave 
(2015), who implemented halo quenching in cosmological hy¬ 
drodynamic simulations by triggering quenching in regions 
dominated by hot (lO^^'A) gas. They reproduced a broad 
range of empirical trends detected in PIO and Woo et al. 
(2013), suggesting that the halo mass remains the determin¬ 
ing factor in the quenching of low-redshift galaxies. 

Another alternative quenching model is the so-called 
“age-matching” prescription of Hearin & Watson (2013) and 
its recently updated version of Hearin et al. (2014). Age¬ 
matching is an extension of the “subhalo abundance match¬ 
ing” (SHAM; Conroy et al. 2006) technique, which assigns 
stellar masses to individual subhalos (including both main 
and subhalos) in the N-body simulations based on halo prop¬ 
erties like the peak circular velocity (Reddick et al. 2013). In 
practice, after assigning M, using SHAM, the age-matching 
method further matches the colours of galaxies at fixed M, 
to the ages of their matched halos, so that older halos host 
redder galaxies. In essence, the age-matching prescription 
effectively assumes a stellar mass quenching, as the colour 
assignment is done at fixed M, regardless of halo mass or en¬ 
vironment, with a secondary quenching via halo formation 
time. Therefore, the age-matching quenching is very simi¬ 
lar to the M.-dominated quenching of PIO, except that the 
second variable is halo formation time rather than galaxy 
environment. 

The key difference between the Mi,- and M,-dominated 
quenching scenarios lies in the way central galaxies become 
quiescent. One relies on the stellar mass while the other on 
the mass of the host halos, producing two very different sets 
of colour-segregated stellar-to-halo relations (SHMRs). At 
fixed halo mass, if stellar mass quenching dominates, the 
red centrals should have a higher average stellar mass than 
the blue centrals; in the halo quenching scenario the two 
coloured populations at fixed halo mass would have similar 
average stellar masses, but there is still a trend for massive 
galaxies to be red because higher mass halos host more mas¬ 
sive galaxies. This difference in SHMRs directly translates to 
two distinctive ways the red and blue galaxies populate the 
underlying dark matter halos according to their M, and M*, 
hence two different spatial distributions of galaxy colours. 

Therefore, by comparing the Wp and AE predicted from 
each quenching model to the measurements from SDSS, we 
expect to robustly distinguish the two quenching scenar¬ 
ios. The iHOD framework we developed in Paper I is ide¬ 
ally suited for this task. The iHOD is a global “halo oc¬ 
cupation distribution” (HOD) model defined on a 2D grid 
of M. and M*, which is crucial to modelling the segrega¬ 
tion of red and blue galaxies in their M, distributions at 
fixed Ml,. The iHOD quenching constraint is fundamentally 
different and ultimately more meaningful compared to ap¬ 
proaches in which colour-segregated populations are treated 
independently (e.g.. Tinker et al. 2013; Rodrfguez-Puebla 
et al. 2015). Our iHOD quenching model automatically ful¬ 
fills the consistency relation which requires that the sum 
of red and blue SHMRs is mathematically identical to the 
overall SHMR. More importantly, the iHOD quenching model 
employs only four additional parameters that are directly 
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related to the average galaxy quenching, while most of the 
traditional approaches require ~20 additional parameters, 
rendering the interpretation of constraints difficult. Fur¬ 
thermore, the iHOD framework allows us to include ~80% 
more galaxies than the traditional HODs and take into ac¬ 
count the incompleteness of stellar mass samples in a self- 
consistent manner. 

This paper is organized as follows. We describe the se¬ 
lection of red and blue samples in Section 2. In Section 3 we 
introduce the parameterisations of the two quenching models 
and derive the iHODs for each colour. We also briefly describe 
the signal measurement and model prediction in Sections 2 
and 3, respectively, but refer readers to Paper I for more 
details. The constraints from both quenching mode analyses 
are presented in Section 4. We perform a thorough model 
comparison using two independent criteria in Section 5 and 
discover that halo quenching model is strongly favored by 
the data. In Section 6 we discuss the physical implications 
of the halo quenching model and compare it to other works 
in 7. We conclude by summarising our key findings in Sec¬ 
tion 8. 

Throughout this paper and Paper I, we assume a ACDM 
cosmology with o'g, h) = (0.26,0.74,0.77,0.72). All the 

length and mass units in this paper are scaled as if the Hub¬ 
ble constant were fOOkm In particular, all the sep¬ 

arations are co-moving distances in units of either /j^'kpc or 
/i^'Mpc, and the stellar mass and halo mass are in units of 
and respectively. Unless otherwise noted, the 

halo mass is defined by Mj, = M 2 oom = 200p„,(47r/3)r2QQ„,, where 
rioom is the corresponding halo radius within which the aver¬ 
age density of the enclosed mass is 200 times the mean mat¬ 
ter density of the Universe, pm- For the sake of simplicity. 
In X- logp X is used for the natural logarithm, and Ig x= logjg x 
is used for the base-fO logarithm. 


2 SAMPLE SELECTION AND SIGNAL 
MEASUREMENT 

In this section we describe the SDSS data used in this pa¬ 
per, especially the selection of the red and blue galaxies 
within the stellar mass samples, and the measurements of 
the galaxy clustering and the g-g lensing signals. We briefly 
describe the overall large-scale structure sample and the sig¬ 
nal measurement, same as that used in Paper I, below in 
Section 2.1 and 2.3, respectively, and refer readers to Paper 
I for details. Here we focus more on the colour cut we em¬ 
ploy to divide the galaxies into red and blue populations in 
Section 2.2. 

2.1 NYU-VAGC and Stellar Mass Samples 

We make use of the final data release of the SDSS (DR7; 
Abazajian et al. 2009), which contains the completed data 
set of the SDSS-I and the SDSS-II. In particular, we ob¬ 
tain the Main Galaxy Sample (MGS) data from the dr72 
large-scale structure sample brightO of the “New York Uni¬ 
versity Value Added Catalogue” (NYU-VAGC), constructed 
as described in Blanton et al. (2005). The brightO sample 
includes galaxies with 10<mr<17.6, where is the r-band 
Petrosian apparent magnitude, corrected for Galactic ex¬ 
tinction. We apply the “nearest-neighbour” scheme to cor¬ 


rect for the 7% galaxies that are without redshift due to 
fibre collision, and use data exclusively within the contigu¬ 
ous area in the North Galactic Cap and regions with angular 
completeness greater than 0.8. The final sample used for the 
galaxy clustering analysis includes 513,150 galaxies over a 
sky area of 6395.49 deg^. A further 5 per cent of the area 
is eliminated for the lensing analysis, due to the absence of 
source galaxies in that area. 

As discussed in Paper I, we further restrict our analysis 
to galaxies above a “mixture limit”, defined as the stellar 
mass threshold above which the galaxy sample is relatively 
complete with a fair mix of red and blue galaxies. The func¬ 
tional form we adopt to describe the mixture limit 
is 

Ig = 5.4 X (2-0.025)“'+ 8.0, (1) 

shown as the thick yellow curves in Fig. 3 (discussed further 
below). By taking into account the sample incompleteness 
in a self-consistent way, the iHOD model is able to model 
the lensing and clustering statistics of all galaxies above the 
mixture limit, ~84% more than the traditional HOD models 
typically include from the same catalogue. 

We employ the stellar mass estimates from the lat¬ 
est MPA/JHU value-added galaxy catalogue*. The stel¬ 
lar masses were estimated based on fits to the SDSS 
photometry following the philosophy of Kauffmann et al. 
(2003) and Salim et al. (2007), and assuming the 
Chabrier (Chabrier 2003) initial mass function (IMF) and 
the Bruzual & Gharlot (2003) SPS model. The MPA/JHU 
stellar mass catalogue is then matched to the NYU-VAGG 
brightO sample. We identify valid, unambiguous MPA/JHU 
stellar mass estimates for all but 32,327 (6.3%) of the MGS 
galaxies. For those unmatched galaxies, we predict their stel¬ 
lar masses using the overall scaling between the two stellar 
mass estimates, depending on the g—r colours (k-corrected 
to z=0.1). 

As sources for the g-g lensing measurement, we use a 
catalogue of background galaxies (Reyes et al. 2012) with a 
number density of 1.2 arcmin^^ with weak lensing shears esti¬ 
mated using the re-Gaussianization method (Hirata & Seljak 
2003) and photometric redshifts from Zurich Extragalactic 
Bayesian Redshift Analyzer (ZEBRA, Feldmann et al. 2006). 
The catalogue was characterised in several papers that de¬ 
scribe the data, and use both the data and simulations to 
estimate systematic errors (see Reyes et al. 2012; Mandel- 
baum et al. 2012; Nakajima et al. 2012; Mandelbaum et al. 
2013). 

2.2 Separating Sample into Red and Blue 

We define quenching by the (^-r)**' colour (after K- 
correction to z=0.1) for three reasons: 1) colour bimodal¬ 
ity is very stable across different environments and red- 
shifts (Baldry et al. 2006); 2) observationally colour is very 
easy to measure robustly, without the need to fit galaxy mor¬ 
phology or brightness profile; and 3) physically colour is the 
result of integrated star formation history, largely immune 
to incidental or minor star formation episodes. In addition, 

' http://home.strw.leidenuiiiv.nl/~jarle/SDSS/ 
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Figure 1. The colour—stellar mass diagram of the SDSS Main Galaxy Sample at four different redshifts (from left to right: z = 0.02, 0.08, 
0.12, and 0.16). The colour maps indicate the co-moving galaxy number density at fixed colour and M*, normalized by the stellar mass 
function at that M^. The red dashed lines are the same across all panels, indicating little redshift evolution in the locus of red sequence 
on the colour-M^ diagram. The black solid lines are the redshift-independent colour cut that we use to divides the galaxies into red and 
blue populations. 
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Figure 2. Illustration of the agreement between using (g-r)® * colour and star formation rate (SFR) as quenching indicators. Each 
panel shows the joint probability density distributions (PDFs) of the colour and the total star formation rate of red and blue galaxies at 
four different stellar masses (from left to right: Ig M,=9.8, 10.2, 10.6, and 11.0). Gray horizontal thick and vertical thin lines indicate the 
colour and the SFR cuts that divide the galaxies into red/blue and passive/active, respectively. The PDFs of red and blue galaxies, each 
normalized to unity, are represented by two separate sets of contour maps with their levels indicated by the corresponding colour bars on 
the right. The red galaxy samples include a very small fraction of dusty star-forming galaxies, with a contamination rate of ~6% - 14% 
depending on stellar mass (marked on the top right of each panel). 


we aim to model and compare to the two separable quench¬ 
ing trends with stellar mass and environment that revealed 
in PIO, who originally chose optical colour as the quenching 
indicator. 

Fig. 1 illustrates the colour-stellar mass dia¬ 
grams (CSMDs) at four different redshifts and the 
stellar mass dependence of the colour cuts we applied 
to divide the red and blue galaxies. In each panel, the 
colour map indicates the distribution of the logarithmic 
comoving number density of galaxies in cells of (g-r)® ' 
and M,, normalized by the stellar mass function (SMF) 
at that M,. This normalization highlights the g-r ranges 
with relatively high concentration of galaxies along the 
colour axis, enhancing the appearance of the “red sequence” 
on each panel. The CSMDs are cut off at different stellar 
masses due to the redshift dependence of the mixture limit, 
which is ultimately related to the flux limit of the spectro¬ 
scopic survey. The red dashed lines going through the red 
sequences are uniform across all redshifts, indicating that 


the loci of the red sequence on the CSMD is independent of 
redshift within our sample. To divide the galaxies into red 
and blue, we therefore define the colour cut to be parallel 
to the red sequence on the CSMD, described by 

= , ( 2 ) 

and indicated by the black solid lines in Fig. 1. The weak 
stellar mass dependence in equation 2 causes a variation of 
(g - between 0.76 and 0.84 within our sample, leading 
to differences in classification between blue vs. red of only a 
few per cent compared to a constant cut at 0.8. 

Whether a galaxy is quiescent or star-forming, how¬ 
ever, is never a clear-cut choice. Galaxy bimodality shows 
in nearly every aspect of galaxy properties, including broad¬ 
band colour, star formation rate (SFR), morphology (e.g., 
late/early-type, De Vaucouleurs/exponential profile), and 
concentration (e.g., Sersic index). Bernardi et al. (2010) 
found that many late-type (Sb and later) galaxies lie above 


MNRAS 000, 1-26 (2015) 






























iHOD quenching model: separating red and blue galaxies 5 



Redshift Redshift 

Figure 3. Selection of the red (left) and the blue (right) galaxy samples on the M,—redshift diagram. In each panel, the color map 
indicates the logarithmic number counts of galaxies at fixed M, and redshift. The thick yellow curve is the “mixture limit” above which 
the completeness of the galaxy samples is relatively high. By selecting almost all galaxies above this mixture limit, our iHOD analysis is 
able to include 84% more than galaxies than the traditional HOD methods. 


the red galaxy colour cut (similar to equation 2) and they 
tend to be edge-on discs reddened by dust. Conversely, some 
early-type galaxies lie below the cut, either showing star¬ 
forming AGN or post-starburst spectrum, with their star 
formation history well described by a recent minor and short 
starburst superimposed on old stellar component (Huang & 
Gu 2009). As a result. Woo et al. (2013) advocated the use of 
SFR as the quenching indicator, and they claimed that one 
third of the red galaxies are star forming. However, the large 
fraction of star-forming contaminants in the “red” popula¬ 
tion in Woo et al. (2013) is mainly caused by the rest-frame 
U—B colour Woo et al. (2013) adopted in defining red galax¬ 
ies, derived from AB magnitudes that are K-corrected from 
the SDSS ugriz photometry. Using the native colour 

largely eliminates the star-formers from the red galaxies. 

Fig. 2 clearly demonstrates the good consistency be¬ 
tween using (g-r)° * colour and SFR as quenching indicators. 
The four panels show the joint 2D probability density dis¬ 
tributions (PDFs) of galaxy colour and the logarithmic SFR 
at four different M, (from left to right: \gM, = 9.8, 10.2, 10.6, 
and 11.0). In each panel, the thick horizontal line represents 
the colour cut defined in equation 2, while the thin vertical 
line indicates the SFR value that saddles the separate SFR 
distributions of passive and active galaxies, which can be 
well described by 

SFR 

Ig = -0.35(lgM. - 10.0) - 10.23, (3) 

in parallel to the star-forming sequence defined in Salim 
et al. (2007). The fraction of dusty, star-forming galaxies 
in the red population decreases from 0.14 to 0.06 as stel¬ 
lar mass increases from lgM,=9.8 to 11.0, significantly lower 
than the one third reported in Woo et al. (2013) using U—B 
colours. Therefore, we conclude that it is robust to use red 
fraction as a proxy for quenching efhciency, and the results 
of our analysis should stay the same if SFR were used for 
the selection of quenched galaxies. 

As described in Paper I, the iHOD model constructs in¬ 
dividual HODs within very narrow redshift slices (we use 


Az=0.01), so that the sample selection does not require a 
single uniform stellar mass range among all the redshift 
slices within that sample, i.e., having a rectangular shape 
on the M,-z diagrams. Fig. 3 illustrates the galaxy samples 
selected on the M,-z diagram within each coloured popula¬ 
tion for the iHOD quenching analysis. The colour intensity 
represents the logarithmic galaxy number counts in cells of 
M, and z- As mentioned in Section 2.1, all selected samples 
have the “wedge”-like stellar mass thresholds described by 
the mixture limit, and thus contain extra galaxies at the far 
end of the redshift range that are usually unused in tradi¬ 
tional HOD analysis. Additionally, since those high redshift 
wedges have a larger comoving volume per unit redshift than 
the low redshifts, they include the most abundant regions on 
the M,-z diagram, corresponding to the reddest regions on 
both panels of Fig. 3. The resultant increase in the selected 
galaxy sample sizes is more than 80% compared to tradi¬ 
tional selections. 

Above the mixture limit, the red galaxies (left panel) 
are two times more abundant than the blue galaxies (right 
panel), despite the fact that the ratio of the two colours is 
close to unity in the spectroscopic survey. Therefore, we can 
afford finer binning in M, in the red galaxy sample than 
in the blue galaxy sample, especially at the high mass end. 
Table 1 summarises the basic information of the two sets 
of sample selections used by the iHOD quenching analysis. In 
total, we divide 228,759 red galaxies and 94,325 blue galaxies 
into eight and six subsamples, respectively. 


2.3 Measuring Galaxy Clustering Wp and 
Galaxy-Galaxy Leasing AE 

We measure the projected correlation function Wp for each 
galaxy sample by integrating the 2D redshift-space correla¬ 
tion function f', 

^\rp, rP)dr^, ( 4 ) 

^ax 
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Table 1 . Red and blue stellar mass bins used for the iHOD quench¬ 
ing analysis, corresponding to the selections in the left and right 
panels of Figure 3, respectively. The red selection includes three 
stellar mass samples above lgM,=ll as opposed to the single 
lgM,>ll sample in the blue selection, while the five lower stel¬ 
lar mass red and blue samples share the same binning in IgM, 
and z. 



Z 

^red 

^blue 

8.5-9.4 

0.01-0.04 

3,224 

10,773 

9.4-9.8 

0.02-0.06 

7,336 

9,356 

9.8-10.2 

0.02-0.09 

28,301 

19,883 

10.2-10.6 

0.02-0.13 

70,514 

29,160 

10.6-11.0 

0.04-0.18 

84,108 

21,058 

11.0-11.2 

0.08-0.22 

22,626 


11.2-11.4 

0.08-0.26 

9,775 


11.4-12.0 

0.08-0.30 

2,875 


11.0-12.0 

0.08-0.30 


4,095 


where r,, and r„ are the projected and the line-of-sight (LOS) 
comoving distances between two galaxies. We measure the 
Wp signal out to a maximum projected distance of = 
20/j^'Mpc, where the galaxy bias is approximately linear. For 
the integration limit, we adopt a maximum LOS distance 
of r““=60/z“'Mpc. We only use the Wp values down to the 
physical distance that corresponds to the fibre radius at the 
maximum redshift of each sample, with fewer data points 
for higher stellar mass (hence larger maximum fibre radius) 
samples. 

The Landy-Szalay estimator (Landy & Szalay 1993) is 
employed for computing the 2D correlation ^‘‘(rp, r„) The er¬ 
ror covariance matrix for each Wp measurement is estimated 
via the jackknife re-sampling technique. We divide the en¬ 
tire footprint into 200 spatially contiguous, roughly equal- 
size patches on the sky and compute the Wp for each of the 
200 jackknife subsamples by leaving out one patch at a time. 
For each stellar mass sample, we adopt the sample mean of 
the 200 subsample measurements as our final estimate of Wp, 
and the sample covariance matrix as an approximate to the 
underlying error covariance. 

For the surface density contrast AE, we measure the 
projected mass density in each radial bin by summing over 
lens-source pairs “Is" and random lens-source pairs “rs", 


Al,(rp) = <Ecrit7f(r,,)> = 


i;,., W/,ef"*Ecrit(Z/,Z.;) 


(5) 


where e, is the tangential ellipticity component of the source 
galaxy with respect to the lens position, the factor of 27? 
converts our definition of ellipticity to the tangential shear 
7 ,, and w/j is the inverse variance weight assigned to each 
lens-source pair (including shot noise and measurement error 
terms in the variance). Ecrit is the so-called critical surface 
mass density, defined as 




4;rG + zi? 

c2 D„ 


( 6 ) 


where A and D, are the angular diameter distances to lens 
and source, and £)/, is the distance between them. We use 
the estimated photometric redshift each source to compute 
Ds and D/j. The factor of (1 + ziY comes from our use of 
comoving coordinates. We subtract off a similar signal mea¬ 
sured around random lenses, to subtract off any coherent 


systematic shear contributions (Mandelbaum et al. 2005); 
this signal is statistically consistent with zero for all scales 
used in this work. Finally, we correct a bias in the signal 
caused by the uncertainties in the photometric redshift us¬ 
ing the method from Nakajima et al. (2012). 

To calculate the error bars, we also used the jackknife 
re-sampling method. As shown in Mandelbaum et al. (2005), 
internal estimators of error bars (in that case, bootstrap 
rather than jackknife) perform consistently with external es¬ 
timators of errorbars for AE on small scales due to its being 
dominated by shape noise. 


3 QUENCHING MODELS AND SIGNAL 
PREDICTIONS 

In this section, we introduce the mathematical descriptions 
of the hybrid (Section 3.1) and the halo (Section 3.2) quench¬ 
ing models. We also briefly describe how to infer the iHODs 
for the red and blue galaxies in Section 3.3, but refer reader 
to Paper I for more details on the iHOD framework. The 
prediction of Wp and AE from each coloured iHOD is rather 
complex but exactly the same as that for the overall galaxy 
populations, therefore we directly refer readers to the rel¬ 
evant sections (4 and 5) in Paper I for details. We ignore 
quenching via mergers in both quenching models consid¬ 
ered below, as merging-induced quenching is negligible at 
z<0.5 (Peng et al. 2010). 


3.1 Hybrid Quenching Model 

The hybrid quenching model parameterizes the red fraction 
as a function of both M, and M*, aiming to mimic the empiri¬ 
cal stellar mass and environment quenching trends observed 
in PIO. In the physical picture implied by this model, ev¬ 
ery quiescent galaxy had spent some portion of its life on 
the star-forming “main sequencing” as a central before the 
eventual quenching (Daddi et al. 2007; Noeske et al. 2007; 
Speagle et al. 2014), due to either the depletion of gas sup¬ 
ply (i.e., stellar mass quenching) or entering another halo as 
a satellite, when environment quenching kicked in. 

While the stellar mass trend is straightforward to pa¬ 
rameterize, it is unclear whether the environment quenching 
trend among satellites can be mimicked by a trend in halo 
mass, as the environment-halo mass relation is very com¬ 
plex and depends strongly on the definition of that environ¬ 
ment. The PIO environment of satellites, as defined by ctsnn, 
shows strong correlation with group richness when the rich¬ 
ness is below five. In richer systems, however, the correlation 
is mostly smeared out and ctj^n instead anti-correlates with 
the halo-centric distance Dg (Peng et al. 2012). This appar¬ 
ent transition between the two richness regimes is caused by 
the increase of D^IRy^^, the ratio between the typical distance 
to the fifth nearest neighbour to the halo virial radius, from 
below to above unity. When Ds > /?vir, cr 5 i,ji,j is roughly pro¬ 
portional to thus more tied to M;,. When Ds < f?vir, 

(Tsnn is essentially an intra-halo overdensity measured at Dg, 
which depends more strongly on Dg than M* due to the steep 
declining slope of the NFW-like halo density profile. 

But for our purposes, what matters is the mean satellite 
quenching efficiency as a function of halo mass Mi ,, averaged 
over galaxies at all Dg within that halo. As pointed out by 
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Figure 4. Illustration of the difference between hybrid (top row) and halo (bottom row) quenching models. For each quenching model, 
we show the red galaxy fraction as a function of both stellar mass and halo mass for the central (left) and satellite (right) populations, 
respectively, predicted from the best-fit parameters. The long arrow in each panel indicates the direction of increasing red fraction, driven 
by either stellar mass (top left), halo mass (bottom left and right), or the combination of both (top right) within the two quenching 
models. 


Woo et al. (2013), the density profile of more massive halos 
falls off less steeply with distance than that of less mas¬ 
sive systems, so the probability of finding the 5-th nearest 
neighbour increases with halo mass. Therefore, the PIO en¬ 
vironment quenching trend can be potentially encapsulated 
within the halo model as a satellite quenching dependence 
on halo mass. 

Assuming stellar mass as the main driver of central 
galaxy quenching, we parameterize the red fraction of cen¬ 
trals as 

,C{M„ M„) = 1 - g{M.) = 1 - exp [- (MJMlf ], (7) 

where Mt is a characteristic stellar mass M/,) = 

(e - l)/e = 0.632) and ji dictates how fast the quenching 
efficiency increases with M,, with /r = 1 being exponential. 
The satellites are subject to an extra halo quenching term 
so that 

M,) = 1 - (8) 

with 

h(M„)^exp[-{Mi,/MlJ], (9) 

where is a characteristic halo mass and v controls the 


pace of satellite quenching. The above equations, including 
g and h as powered exponential functions, are very similar 
to the fitting formula adopted in Baldry et al. (2006) and 
Peng et al. (2012). 

The top left and right panels of Fig. 4 illustrate the 
central and satellite red fractions, computed from the best-fit 
hybrid quenching model via Equations 7 and 8, respectively. 
The arrow in each panel points in the direction of increasing 
quenching efficiency on the 2D plane of M, and Mi,. For 
the central galaxies, although the quenching is driven by 
M, along the horizontal axis, the red fraction still shows 
strong increasing trend with Mi, due to the tight correlation 
between M, and Mi,, i.e., the SHMR of central galaxies. The 
2D distribution of satellite red fractions displays a “boxy” 
pattern, echoing the separate stellar mass and environment 
quenching trends detected in PIO. 

Combining the central and the satellite terms, the red 
fraction of galaxies with stellar mass M, inside halos of total 
mass Ml, is 

r\M„ Ml,) = fUM„ Ml,) fl^(M„ Ml,) + 

[1 - /.a,(A7„ M„)] f^(M„Mi,), (10) 

where fsM(M,,Mi,) is the satellite fraction that can be pre- 
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dieted by the overall iHOD model from Paper I. For the hy¬ 
brid model, equation 10 can be reduced to 

( 11 ) 


3.2 The Halo Quenching Model 

As described in the introduction, the halo quenching model 
relies on halo mass alone to quench both central and satel¬ 
lite galaxies, and Gabor & Dave (2015) demonstrated that 
it also naturally explains the stellar mass and environment 
quenching trends seen in PIO, by embedding galaxies in 
massive halos filled with hot gas via virial heating. How¬ 
ever, depending on the exact physical processes driven by 
M;,, halo quenching may apply to the central and satellites 
differently. For instance, while the central galaxies in halos 
above Mshock could be quenched by shocked-heated gas and 
then maintain a high gas temperature via the “raido”-mode 
feedback from AGNs, the satellite galaxies in the those ha¬ 
los may still retain some cold gas as the “central” galaxies of 
its own coherent sub-group. Therefore, the satellite galaxies 
continue to accrete gas and convert it to stars over a pe¬ 
riod of ~1 Gyr after entering into a larger halo (Simha et al. 
2009). Similar processes like slow strangulation (assuming no 
accretion onto satellites) also produce prolonged quenching 
actions (Peng et al. 2015). In this case, the halo quenching 
of centrals and satellites are somewhat decoupled, and the 
quenching of satellites is a more gradual process than that 
of centrals. 

Therefore, unlike the hybrid model, we describe the red 
fractions of centrals and satellites as two independent func¬ 
tions of Ml,-. 


//:/(44„44„)= l-/“r(44/.) 


= l-exp[-(44*/44f)''‘]. 

(12) 

and 


fllt{M„Mi,)=l-f^\Mi,) 


= l-exp[-(44„/44ff]. 

(13) 

where Mf and 44^“ are the critical halo 

masses responsi- 


ble for triggering quenching of central and satellites, respec¬ 
tively, and and are the respective powered-exponential 
indices controlling the transitional behavior of halo quench¬ 
ing across the critical halo masses. 

Similarly, the bottom left and right panels of Fig. 4 illus¬ 
trate the central and satellite red fractions, computed from 
the best-fit halo quenching model via Equations 12 and 13, 
respectively, with arrows indicating halo mass as the sole 
driver for quenching in both populations. The orthogonal¬ 
ity of the hybrid and halo quenching directions for central 
galaxies is the key distinction that we look to exploit in this 
paper, by identifying its imprint on the clnstering and g-g 
leasing signals of red and blue galaxies. The total red frac¬ 
tion can be obtained by substituting equations 12 and 13 
into equation 10. 

Finally, we emphasize that in reality the true quenching 
arrow could be pointing anywhere between the two orthog¬ 
onal directions, i.e., a more generalized quenching model 
consisting of a linear mixture of the two, with the linear 


coefficients varying as functions of M, and Mi, as well — 
schematically, 

Q\mt — 44/,) X ^hybrid "t (1 ~ n)(A4*, Mj,')) X 0halo- (1^) 

However, as a first step of constraining quenching, the goal 
of this paper is to find out if u) is closer to zero (i.e., halo¬ 
quenching dominated) or unity (i.e, hybrid-quenching dom¬ 
inated) . 


3.3 Prom Quenching Models to Colour-segregated 
iHODs 


In order to predict the iVp and AE for the red and blue galax¬ 
ies in each quenching model, we construct iHODs for both 
coloured populations by combining the overall iHOD with 
f'^'^iM,, Ml,) predicted by that quenching model. 

Let us start with the red galaxies. The key is to derive 
p^’^'^iM,, Ml,), the 2D joint PDF of the red galaxies of stellar 
mass M, sitting in halos of mass Mi,, given the 2D PDF of 
the overall galaxy population p(M,,Mi,) inferred from Paper 

I, 

/'“(M,, Ml,) = (15) 

/tot 


where is the overall red faction of all galaxies, obtained 
via 



f^'^{M„Mi,)p{M„ Ml,) dMi,dM,. 


(16) 


As described in Paper I, iHOD predicts the Wp and AE 
signals for a given galaxy sample by combining the predicted 
signals from individual narrow redshift slices, each of which 
is described by a single standard HOD. 

For deriving standard HODs within redshift slices for 
the red galaxies, we need 




M,) 

p’^’^'^iM,) 


(17) 


while is the predicted parent (i.e., including observed 

and unobserved galaxies) SMF of the red galaxies normal¬ 
ized by their total number density 






n 


red 

tot 





(18) 


Finally, we arrive at the HOD of red galaxies at any 
redshift z as 






p"‘‘(44/,|44,)®:,t(A4,k)d44„ 


(19) 


where <l)[)j(*^(44.|z) is the observed SMF of red galaxies at red¬ 
shift z, directly accessible from the survey. For modelling the 
samples defined in Figure 3 for the iHOD analysis, we mea¬ 
sure the observed galaxy SMF at each redshift, and then 
obtain the HOD for that redshift slice using Equation (19). 
In this way, we avoid the need to explicitly model the sample 
incompleteness as a function of 44, and/or 44/,. 

For the blue galaxies, we apply the same procedures 
above to obtain ^A''’’“‘^(44;,|z)^ from p'’'“‘=(44/,, 44,), by substitut¬ 
ing /‘■'=‘‘(44., 44,.) with /““(A4„ 44,,)=1 - /‘'"‘‘(Af,, 44*) in equa¬ 
tions 15-19. 

Fig. 5 illustrates the two sets of coloured iHODs derived 
from the best-fit hybrid (top row) and halo (bottom row) 
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Figure 5. The iHOD, i.e., the average number of galaxies per dex in stellar mass within halos at fixed mass of the red (left column) 
and blue (right column) populations, predicted from the best-fit hybrid (top row) and halo (bottom row) quenching models based on 
the same iHOD of the overall population. Black and white contour lines in each panel indicate the iHOD distributions for the central 
and satellite galaxies, respectively, while the colour contours show the iHOD for the total galaxy population. The two best-fit models 
produce similar overall patterns of galaxy occupation, with very subtle but important differences — compared to hybrid quenching, halo 
quenching produces stronger segregation between red and blue centrals in high vs. low mass halos. 


quenching models. In each row, the left and right panels dis¬ 
play lg(dAi(M,|M;,)/dlgM,), the average log-number of galax¬ 
ies per dex in stellar mass within halos at fixed mass, for 
the red and blue populations, respectively. The white and 
black contour lines highlight the central and satellite galaxy 
occupations separately on the plane. All panels re¬ 

veal the same generic pattern, consisting of a tight sequence 
that corresponds to the SHMR of the central galaxies, and 
a cloud underneath occupied by the satellite galaxies. The 
level of similarity exhibited by the red galaxies is especially 
high between the two quenching models (left column). 

However, comparing the left and right panels in the 
same row (i.e., red vs. blue galaxies in the same quench¬ 
ing model), the red centrals are more preferentially sitting 
in the high-M, and high-M;, region than in the low-M, and 
low-M;, one, while the opposite is true for the blue centrals. 
This segregation happens regardless of quenching models, 
confirming our notion that it is difficult to unambiguously 
disentangle the two quenching directions, despite their or¬ 
thogonality, by merely examining the quenching trend with 
M,, or some surrogate of M;, that has substantial scatter 
about the true M/, (e.g., group richness). 


The satellites are quenched by M;, in both models, but 
are also partially by M, in the hybrid model. Thus, there are 
more high-M, blue satellite galaxies within massive halos in 
the halo quenching model (bottom right panel) than in the 
hybrid model (top right panel). In addition, the low mass 
halos in the hybrid quenching model are more likely to host 
blue dwarf satellites than in the halo quenching model. 


The two sets of iHOD models, presented in this Section 
and in Fig. 5, are the analytical foundation that allow us to 
predict the Wp and AE signals as functions of the four pa¬ 
rameters in each quenching model. Any difference shown in 
Fig. 5 between the two models will be propagated to the dif¬ 
ferent behaviours in the final predictions of Wp and AX, and 
is thus detectable by comparing the two sets of predictions 
to the measurements from SDSS galaxies. 
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Table 2. Description, prior specifications, and posterior constraints of the parameters in the halo (top) and hyrbid (bottom) quenching 
models. All the priors are uniform distributions running across the entire range of possible values for the parameters, and the uncertainties 
are the 68% confidence regions derived from the ID posterior probability distributions. 


Parameter 

Description 

Uniform Prior Range 

prior case 

fixed case 

Halo Quenching Model 

\gMf [ft-'Mo] 

IgMf 

Characteristic halo mass for central galaxy quenching 

Pace of central galaxy quenching with halo mass 
Characteristic halo mass for satellite galaxy quenching 

Pace of satellite galaxy quenching with halo mass 

[11.0, 15.5] 

[0.0, 3.0] 

[11.0, 15.5] 

[0.0, 3.0] 

12 

i/,./,u_oo8 
n ”10+0.04 

9-JO_003 

12 17+*'32 

^•^'^-0.02 

12 

iZ.zj_oo6 

0 42+0 03 
^■^^-0,03 

0 16+0-03 
-0.03 

Hybrid Quenching Model Qhybrid 

\gMi [r^Mol 

lgM« [r‘Mo] 

V 

Characteristic stellar mass for central and satellite quenching 
Pace of galaxy quenching with stellar mass 
Characteristic halo mass for satellite galaxy quenching 

Pace of satellite galaxy quenching with halo mass 

[9.0, 12.0] 

[0.0, 3.0] 

[11.0, 15.5] 

[0.0, 3.0] 

10 50'^^ ^ 

-0,06 

13 76+^ 

-0,14 

0 15+0.05 

10 55+0 03 

^■^^-0,05 

13.63)^ 


u 
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Figure 6. Confidence regions from our halo-quenching model 
analysis of the clustering and lensing of the red and blue galaxies 
in the 2D planes that comprised of all the pair sets of the four 
quenching parameters. Histograms in the diagonal panels show ID 
posterior distributions of individual parameters. Contour levels 
run through confidence limits of 95 per cent (light brown) and 
68 per cent (dark brown) inwards. The filled contours show the 
constraints from our fiducial model where the iHOD parameters 
that describe the overall galaxy population are drawn from priors 
informed by the analysis in Paper I, while the open contours in 
the panels of the upper triangle are the constraints from a simpler 
model where the overall iHOD parameters are kept fixed at the 
best-fit values derived from Paper 1. 

4 CONSTRAINTS ON THE TWO 
QUENCHING MODELS 

4.1 Constraints of the Quenching Parameters 

Ideally one would constrain both the iHOD parameters and 
the quenching parameters together, by simultaneously fit¬ 
ting to the Wp and AS measurements of the overall, red, 
and blue galaxies. However, since the measurements of the 
overall population have the highest signal-to-noise ratio and 
the overall iHOD does not include quenching, it is concep- 
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Figure 7. Similar to Fig. 6, but for the hybrid quenching model. 


tually more reasonable to adopt a two-step scheme. In the 
hrst step we constrain the iHOD parameters using only the 
measurements of the galaxy samples without dividing by 
colour (i.e., Paper I). In the second step, when constraining 
the quenching parameters, we either fix the best-fit iHOD pa¬ 
rameters (i.e., the “fixed case”) or input the iHOD constraints 
from Paper I as priors (i.e., the “prior case”). In particular, 
for the prior case we draw the iHOD parameters from the 
joint prior distribution represented by the Markov Chain 
Monte Carlo (MCMC) samples derived in Paper I. For each 
quenching model, we adopt the results from the prior case 
as our fiducial constraint in the following analysis. 

In addition to the powerful statistical features of the 
iHOD framework inherited from Paper I, our quenching anal¬ 
ysis also adds two important advantages compared to the 
traditional HOD modelling of red and blue galaxies. Firstly, 
traditional HOD studies of red and blue galaxies treat the 
two populations independently, so that the total number of 
HOD parameters inevitably doubles compared to the mod¬ 
elling of the overall galaxy population (e.g.. Tinker et al. 
2013; Rodrfguez-Puebla et al. 2015). In our analysis, the 
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Figure 8. Comparison between the measurements of clustering (top row) and lensing (bottom row) signals to the predictions by the two 
best-fit quenching models, for the eight red stellar mass samples (left column) and the six blue samples (right column). In each panel, the 
signals from bottom to top are in ascending order of the average stellar mass of the galaxy samples, each shifted from its adjacent sample 
by 1.0 dex to avoid clutter. Points with errorbars are the measurements, while the solid and the dashed curves are the predictions from 
the best-fit halo and hybrid quenching models, respectively. The two best-fit models provide equally adequate fits to the measurements 
for the red galaxies, but predict different g-g lensing signals of high-M, blue galaxies (bottom right; the two high mass samples). 


red and blue populations are derived not from scratch, but 
by filtering the overall iHOD with the red fraction predicted 
by each quenching model, which is described by only four 
simple yet physically meaningful parameters. Our method 
also guarantees that the sum of the red and blue SHMRs 
is mathematically identical to the overall SHMR. Secondly, 
the traditional method usually parameterizes the red galaxy 
fraction as a ID function of halo mass, while our method 
affords a 2D function of defined on the M,-Mi, plane, 
which is crucial to the task of examining stellar mass as a 
potential driver for quenching. 

For each quenching model, we infer the posterior prob¬ 
ability distributions of the four model parameters from the 
Wp and AE measurements of the eight red and six blue galaxy 
samples within a Bayesian framework. We model the combi¬ 
natorial vector X of the Wp and the AE components of the red 
and blue galaxies as a multivariate Gaussian, which is fully 


specified by its mean vector (x) and covariance matrix (C). 
The Gaussian likelihood is thus 


£{x\e) = |Cr‘/" exp 


(x-x)^C-‘(x-x)\ 

2 j’ 


( 20 ) 


where 


e = {lgMt,fi,lgMl,v] (21) 

in hybrid quenching, and 

e = {lgMf,/,lgMf,;/^l (22) 

in halo quenching. 

We adopt flat priors on the model parameters, with a 
uniform distribution over a broad interval that covers the 
entire possible range of each parameter (see the 2nd column 
of Table 2). The final covariance matrix C is assembled by 
aligning the error matrices of Wp and AE measured for in¬ 
dividual coloured samples along the diagonal blocks of the 
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Figure 9. Comparison between the clustering (left) and lensing (right) signals predicted by the halo (solid) and hybrid (dashed) 
quenching models for the red (red), and blue (blue) galaxies. The thick gray curves indicate the best-fit iHOD model prediction for the 
overall galaxy population, with the stellar mass ranges marked on the right panel. The two quenching models generally predict similar 
scale and stellar mass dependences of the red and blue deviations from predictions of the overall population (comparing solid and dashed 
lines of the same colour), except for the high mass blue lenses (comparing the blue solid lines to the blue dashed lines in the right panel). 


full N X N matrix. We ignore the weak covariance between 
Wp and AE (with the covariance being weak due to the fact 
that AE is dominated by shape noise), and between any two 
measnrements of the same type but for different stellar mass 
or coloured samples. 

Fig. 6 presents a summary of the inferences from the 
halo quenching model analysis, showing the ID posterior 
distribution for each of the four model parameters (diag¬ 
onal panels), and the 95% and 68% confidence regions for 
all the parameter pairs (off-diagonal panels). In the pan¬ 
els of the lower triangle, we highlight the results from our 
fiducial model, i.e., the prior case, employing the iHOD pa¬ 
rameter constraints from Paper I as priors. In each panel 
of the upper triangle, we compare the constraints from the 
fiducial analysis (filled contours) to that of the fixed case 
analysis, which keeps the iHOD parameters at their best-fit 
values derived from Paper I. The two analyses are consis¬ 
tent with each other, implying that the explanation of the 
red and blue signals does not require any modification in 
the description of the overall galaxy population. The two 
inferred characteristic halo mass scales are very similar to 
the critical shock heating mass scale, while 

the two powered-exponential indices, and yu‘, indicate that 
the central and the satellite quenching transition differently 
across that shared characteristic halo mass. We defer the 
detailed discussion of the physical implications of the halo 
quenching constraints in Section 6. The 68% confidence re¬ 
gions of the ID posterior constraints are listed in Table 2. 

Similarly, Fig. 7 presents the constraints on the hybrid 
quenching model. The critical stellar mass for quenching all 
galaxies is M’=3.16(±0.31) x echoing the char¬ 

acteristic stellar mass for downsizing at the low redshift. 
The stellar mass quenching index p is slightly below unity, 
the value required for maintaining the observed redshift- 


independence of Schechter M" and faint-end slope of the 
star-forming galaxies in the stellar mass quenching formal¬ 
ism proposed in PIO. The characteristic halo mass for the 
quenching of satellites is much higher than albeit with 
a similar quenching index of v=0.15 ± 0.05. 


4.2 Best-fit Model Predictions 

Fig. 8 compares the clustering (top row) and g-g lens¬ 
ing (bottom row) signals measured from SDSS (points with 
errorbars) to those predicted by the best-fit halo (solid lines) 
and hybrid (dashed line) quenching models, for the eight 
red (left column) and the six blue (right column) stellar mass 
samples. In terms of the overall goodness-of-fit, the best-fit 
halo quenching model yields a. of 701.0, while the hybrid 
quenching model has a worseof 736.8. The reducedval¬ 
ues are thus 1.60 and 1.68 for the halo and hybrid quenching, 
respectively, both providing reasonable fits to the data, con¬ 
sidering that the uncertainties in the measurements of the 
low-M, samples are under-estimated. We defer a discussion 
of the statistical significance of both best fits to the upcom¬ 
ing section. 

For the red galaxy samples, the two quenching models 
predict very similar signals except for the two lowest stel¬ 
lar mass bins. Unfortunately the Wp measurements in these 
two bins are severely affected by the underestimated cosmic 
variance due to the small volumes, with highly correlated un¬ 
certainties on all scales. Therefore, neither quenching model 
gives an adequate fit to their Wp signals. The AE signals of 
the two lowest mass bins are less affected by cosmic variance 
because the measurement error is dominated by shape noise, 
and are thus better described by the two quenching model 
predictions. 

The two quenching models also predict very similar sig- 
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Figure 10. Decomposition of the clustering (top) and lensing (bottom) signals of the high mass blue galaxy sample (Ig M,=[11.0,12.0]), 
by the best-fit hybrid (left) and halo (right) quenching models. In each panel, the green curve shows the predicted signal for the overall 
galaxy population at this stellar mass. The measurement and the predicted signal for blue galaxies are shown by the points with errorbars 
and the thick blue curve, respectively, with the reduced marked on the top. The predicted signal is then decomposed into different 
components involving separate contributions from centrals and satellites (see text for details). In particular, the “l-halo non-satellite” 
terms in the AE panels indicate the underlying dark matter density profiles of the host (sub)halos. The halo quenching model provides 
a much better fit to the Wp and AE measurements by allowing substantially fewer massive halos to host these high-M, blue galaxies as 
centrals, and a higher satellite fraction among these massive blue galaxies. 


nals for the blue galaxies, except for the high mass ones 
with lgM,>ll. While both quenching models give adequate 
fits to the Wp signals of these massive blue galaxies, the halo 
quenching model produces much better fit to their AE sig¬ 
nals than the hybrid quenching model, driving most of the 
difference in the log-likelihoods (i.e., the values) of the 
two best-fit models. This difference revealed by the massive 
blue galaxies, as will be discussed further later, is the key to 
distinguishing the two quenching models. 

Fig. 9 highlights the split between the red and blue 
galaxies from the overall population in the Wp (left) and 
AE (right) signals, predicted by the two best-fit quenching 
models for the eight stellar mass bins marked in the right 
panel. In each panel, the thick gray curves are the iHOD pre¬ 
dictions for the overall galaxy samples, which bifurcate into 
the thin red and blue curves, i.e., predictions for the red 
and blue galaxies. Solid and dashed line styles indicate the 
halo and hybrid quenching models, respectively. As seen in 
Fig. 8, the two quenching models predict very similar bi¬ 
furcation signatures, except for the high mass bins where 
the hybrid quenching predicts a stronger large-scale bias, a 
weaker small-scale clustering strength, but a stronger small- 
scale g-g lensing amplitude, than the halo quenching for the 
blue galaxies. Unfortunately the measured Wp signals for the 


high mass galaxies are cut off at small scales due to fibre col¬ 
lision, and the measurement uncertainties in the large-scale 
Wp are not small enough to distinguish the two quenching 
predictions (top right panel of Fig. 8). 

Therefore, the g-g lensing of the massive blue galax¬ 
ies clearly provides the most discriminative information, as 
shown in the right panel of Fig. 9. For blue galaxies above 
the halo quenching model predicts substantially 
lower weak lensing amplitudes than the hybrid model on all 
distance scales, and thus provides a much better fit to the 
measurements (see bottom right panel of Fig. 8). 

To understand the discrepancy between the two quench¬ 
ing predictions for the massive blue galaxies, we show the 
decomposition of Wp (top row) and AE (bottom row) signals 
predicted by the best-fit hybrid (left column) and halo (right 
column) quenching models for the lgM, = 11.0-12.0 blue sam¬ 
ple in Fig. 10. In each panel, the data points with error- 
bars and the thick blue curve are the measured and pre¬ 
dicted signals for the blue sample, while the iHOD predic¬ 
tion for the overall sample is shown by the thin green curve. 
The best-fit quenching model prediction is then decomposed 
into contributions from the 1-halo and 2-halo (thin dot¬ 
ted) terms. For Wp the 1-halo term includes the contribu¬ 
tions from centrel-satellite pairs (thin solid; “1-h c-s”) and 


MNRAS 000, 1-26 (2015) 




















14 Zu & Mandelbaum 2015 


satellite-satellite (thin dashed; “1-h s-s”) pairs; For AE the 
1-halo term consists of a satellite term (thin dashed) and 
a non-satellite term (thin solid). We also include a point 
source stellar mass term in AE, which is model-independent 
and negligible on most of the relevant scales (not shown 
here). Most importantly, the 1-halo non-satellite term is di¬ 
rectly related to the average dark matter density profile of 
the host halos (including both the main halos for centrals 
and the subhalos for satellites), and its amplitude is propor¬ 
tional to the average mass of those halos. The halo quenching 
model clearly provides a much better £t to the data than the 
hybrid model, with factors of two and four improvement in 
for Wp and AE, respectively. In addition, Fig. 10 shows the 
crucial advantage of including g-g lensing in the joint analy¬ 
sis — since the best-fit hybrid quenching model adequately 
describes the galaxy clustering (x^!N .Ql), its deficiency 
would not be exposed unless we compare the AE predictions 
to data IN='i.%A). 

Fig. 10 reveals two major differences between the two 
quenching model predictions: 1) the halo quenching model 
predicts a much higher satellite fraction among the massive 
blue galaxies than the hybrid model, hence the more promi¬ 
nent “1-halo satellite” terms; and 2) the average (sub)halo 
mass of those massive blue galaxies predicted by the halo 
quenching model is much lower compared to the hybrid 
model prediction, hence the lower g-g lensing amplitudes 
and better fit to the data. Roughly speaking, since the hy¬ 
brid quenching model relies on the stellar mass to quench 
central galaxies, it tends to place central galaxies at fixed 
M, into similar halos regardless of their colours. However, in 
the halo quenching model any galaxies that are unquenched 
have to live in lower mass halos than their quenched coun¬ 
terparts with similar M,. In the section below we will argue 
that, the discrepancy between the average halo masses of 
the massive blue galaxies predicted by the two quenching 
models is insensitive to the details in the model parameters, 
therefore can be used as a robust feature for identifying the 
dominant quenching driver. 


4.3 Origin of Host Halo Mass Segregation 
between Red and Blue Centrals 

Comparison between the two best-fit predictions in Sec¬ 
tion 4.2 reveals that {Mi,\M,), the average host halo mass 
at fixed stellar mass (i.e., the mean halo-to-stellar mass re¬ 
lation; HSMR), is potentially the key discriminator of the 
two types of quenching models. In particular, by predicting 
a lower for the blue centrals, the halo quenching 

model provides a much better fit to the Wp and AE signals of 
the massive blue galaxies than the hybrid quenching model. 
But before going any further, we need to understand the 
cause of this discrepancy between the two quenching models, 
especially to answer the following questions. Firstly, what is 
the origin of the host halo mass segregation between the two 
colours? Secondly, is the halo quenching necessary for pre¬ 
dicting the strong segregation in between the red 

and blue centrals, and can the stellar mass quenching pro¬ 
cess produce an equally low halo mass for the massive blue 
centrals with a different /r? 

For red or blue central galaxies, the conversion from 
the mean SHMR (i.e., <M,|M;,)) to its inverse relation, the 
HSMR, is highly non-trivial. Using the blue centrals as an 


example, the HSMR can be computed from 


where 






(23) 


(24) 


In the above equation, is the PDF of blue cen¬ 

tral galaxy stellar mass at fixed M/,, determined by the mean 
SHMR of the blue centrals and its scatter, /““ is the blue 
fraction of centrals, and du/dM* is the halo mass function. 
Therefore, for given cosmology the HSMR of the blue central 
galaxies has two components, the blue central SHMR (both 
mean and scatter) and the blue fraction of centrals. To un¬ 
derstand for both colours more quantitatively, we 

start by examining the red and blue SHMRs predicted by 
the two models. 

The top and bottom panels in the left column of Fig. 11 
show the mean SHMRs of the total, red, and blue central 
galaxies, predicted by the best-fit hybrid and halo quenching 
models, respectively. Coloured bands indicate the logarith¬ 
mic scatters about the mean relations. The hybrid quenching 
model predicts a segregation in M, between the red and blue 
central galaxies at fixed halo mass, as the high M, galaxies 
are more likely to be quenched. The halo quenching model, 
however, predicts exactly the same SHMRs for all three pop¬ 
ulations, as galaxies at hxed halo mass are equally likely to 
be quenched regardless of stellar mass. The red and blue 
segregation in M,, or the lack thereof, is best illustrated in 
the two right panels of Fig. 11, using three halo masses as 
examples (lg4/, = 11.5, 12.5, 14.0). 

In each panel, the total filled area for each halo shows 
the stellar mass distribution of central galaxies in that halo. 
The width of the distribution decreases with halo mass due 
to the flattening of SHMR on the high mass end. Under 
each distribution, the red and blue shaded areas represent 
the contributions from the red and blue centrals, so that the 
sum of the red and blue SHMRs exactly recovers the total 
SHMR. In the hybrid quenching model for any given halo 
mass, the red galaxy distributions are shifted to higher M, 
compared to the blue distributions, which are indicated by 
the dashed histograms and are equivalent to the blue shaded 
regions. The halo quenching model produces zero such shift. 
The non-zero shift in the hybrid model drives the SHMR 
of the blue central galaxies to become shallower than that 
of the red centrals as seen in the top left panel of Fig. 11. 
Naively, one might think that this shift will also cause the 
blue centrals to reside in more massive halos than the red 
centrals if we simply compare the inverse functions of the 
two SHMRs — a shallower (blue) SHMR maps the same M, 
on the y-axis to a higher halo mass on the x-axis. 

However, a more careful inspection of the segregation 
patterns reveals a second, and much more important differ¬ 
ence in the predicted fraction of blue galaxies among cen¬ 
trals — blue centrals persist in all halo masses in the hybrid 
quenching model, but barely show up in the ha¬ 

los in the halo quenching model. The left panel of Fig. 12 
illustrates the blue fractions as functions of M;, predicted by 
the best-fit halo (thick black solid) and hybrid (thick blue 
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Figure 11. Comparison between the stellar-to-halo mass relations (SHMRs) derived by the hybrid (top) and halo quenching (bottom) 
models. On the left panel of each row, red solid, blue solid, and green dashed curves indicate the average logarithmic stellar mass of red, 
blue, and any central galaxies at fixed halo mass, respectively, predicted by the corresponding quenching model. The shaded band of 
each color shows the log-normal scatter about the average relation. The three colored ticks on the top indicates the three halo masses for 
which we show the central galaxy stellar mass distributions in the right panel. The integral of the distribution for each halo mass (shown 
in the legend) is the expected total number of central galaxies at that halo mass, which can be decomposed into the contributions from 
red (red shaded histogram) and blue (blue dashed histogram or blue shaded area) centrals. The halo quenching model produces exactly 
the same SHMRs for the two coloured populations as the overall relation, due to the lack of correlation between quenching and stellar 
mass at fixed halo mass, whereas the hybrid quenching model makes fractionally more red centrals at higher M,. 


dashed) quenching models. The amplitude of in hybrid 
quenching also depends on M, and the blue dashed curve is 
the average blue fraction over all galaxies above lO^h^^MQ. 
While in the halo quenching case strictly follows the 
powered exponential form (i.e., equation 12), in the hybrid 
case it is affected by both the stellar mass quenching and 
the slope of the SHMR. We also show the blue fraction of 
satellites (thin solid) derived by the halo quenching analysis, 
which exhibits a slower decline with Mi, compared to that of 
centrals. 


The blue galaxy fractions of centrals decline rapidly 
with halo mass in both quenching models, but the speed 
of decline varies differently as a function of halo mass be¬ 
tween the two models. To investigate this quantitatively, we 
define the “central galaxy quenching rate” as a function of 
halo mass, as the logarithmic rate at which the blue 

fraction declines with halo mass, d In/;!),'“/d In M;,, which is 
shown in the right panel of Fig. 12 for each model. As ex¬ 
pected, the halo quenching produces a steady increase of \e\ 


with Ml, 


^halo — 


/ dlnC° \ 
I dlnM* 


M'i 


(25) 


For the hybrid quenching case, f^en^iM;,) experiences a rapid 
decline at low Mi, and then a gradual one at high Mi,. This 
shift in gear can be understood as follows. The central galaxy 
quenching rate ehybrid depends on both the f^^‘(Mf) and the 
derivative of the SHMR, so that 


/dl^\ _ / dln^»r \ / dlnM. \ 

I d In Ml, 1, d In M. i d In mJ ' 


(26) 


Since /ce*n°(lH,) also has a powered exponential form (see 
equation 7), 


/ dln/“r \ 

I dlnM, 


OC 


(27) 


The slope of the SHMR dlnM./dlnM;, is tightly con¬ 
strained by Paper I, which described the SHMR 
/sHMR= exp (In M.(Ma)) as the inverse of 

F 


In — = B In m -I- 
Mi ^ 


m 


I + m 


(28) 
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Figure 12. Left panel: blue galaxy fractions as functions of halo mass predicted by the two best-fit quenching models. Black and Blue 
thick solid curves are the blue fractions within centrals /ce*n°(M;i) predicted by the best-fit halo and hybrid quenching models, respectively. 
Black thin curve indicates the blue satellite fraction predicted by the best-fit halo quenching model. Right panel: the magnitude of the 
logarithmic slopes of characterised by \e{M),)\ in the two quenching models. The hybrid quenching model predicts a rapid decline 

of at low masses and then slows drastically at high masses, while the halo quenching model maintains a steady decline over all 
masses. The two dotted lines are the asymptotic behaviors of the hybrid model at very low and high masses, derived in equation 31. The 
slow decline of /cen° in the hybrid model over-predicts the fraction of blue centrals among massive halos. 


where m = x and M{~\3 x 

are the characteristic stellar and halo mass that 
separate the behaviours in the low and high mass ends, and 
the remaining parameters control the running slopes of the 
SHMR. Assuming reasonable values of the slope parame¬ 
ters (i.e., /3~0.33, S~0A2, 7~1.21; see Paper I), equation 28 
can be approximated by 


j Ml, j /31nra-0.5, m«l 

^ Ml \ m''-I-IjSlnm - 0.5], m»l. 


(29) 


Clearly, the SHMR is a steep power-law relation at the low 
mass end, with M, oc ~ M^, whereas at the high mass 
end the slope of SHMR is very shallow, with M, oc (In Mi,yA ^ 

Therefore, the slope of the SHMR is 


d In M, ( 1 IP, Ml, « Ml 

dlnM,, ~ I ((51nM;.)^‘, M;. » Mi. 


Combining equations 26, 27, and 30, we arrive at 


fhybrid 


-Mf, 


Ml, « Ml 
Ml, » Ml. 


(31) 


halos as efficiently as the hybrid model, but by maintaining 
a steady quenching rate at e{Mi,) -the halo quench¬ 

ing model is able to quench almost all centrals in the very 
massive halos. The hybrid quenching model, however, is rel¬ 
atively inefficient to quench massive central galaxies in the 
very high mass halos. When calculating the HSMR using 
equation 24, this difference in e completely dominates the 
effect due to the slight difference between the two coloured 
SHMRs. Therefore, the stellar mass quenching, due to its 
slow central galaxy quenching rate on the high mass end, is 
incapable of producing a strong segregation in the HSMR 
between the two colours. In order for the hybrid quench¬ 
ing model to mimic the steeper slope of ehaio(M;,), the stellar 
mass quenching trend would have to drop so precipitously 
that the abundance of blue galaxies is cut off beyond some 
maximum stellar mass, which is ruled out by the observed 
SMFs of blue galaxies (see Fig. 14). Therefore, we further 
emphasize that this slow quenching rate with halo mass in 
the hybrid model is caused by the changing slope of SHMR 
across Mi, and is thus insensitive to the stellar mass quench¬ 
ing prescriptions, e.g., the value of p. 


Assuming p = 0.67 from the best-fit hybrid quenching model, 
we have 


( -Ml, Ml, « M| 

-(lnM/.)“6, Ma»Mi. 


(32) 


The above equation is shown as the dotted lines on the right 
panel of Fig. 12, roughly reproducing the two distinctive 
asymptotic behaviours of ei,ybrid at the low and high mass 
ends. The actual slope of ehybrid is steeper than predicted by 
equation 32 at high masses, where equation 30 becomes less 
accurate. 

The comparison between ehaio and ehybrid in Fig. 12 (i.e., 
equations 32 and 25) clearly reveals that, the halo quench¬ 
ing model does not quench central galaxies in the low mass 


To summarize the findings above using the quenching 
diagram of Fig. 4, the steep slope (M,~Mj) of the SHMR 
below Ml makes the SHMR more aligned with the quench¬ 
ing arrow along the M,-axis (i.e., stellar mass quenching, see 
top left panel of Fig. 4), causing progressively more galaxies 
to be quenched at higher halo mass. Above Mi, however, 
the SHMR becomes shallower {M,~(ln Mi,y *) and is almost 
perpendicular to the quenching arrow, leaving a substantial 
number of blue centrals in massive halos. As a result, the 
massive blue centrals are extremely scarce in the 
halos in the halo quenching model (bottom right panel of 
Fig. 11), but have a much stronger presence within those 
halos in the hybrid quenching model (top right panel of 
Fig. 11). By the same token, a strong segregation in the host 
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halo mass between the red and blue centrals would points 
to the necessity of a dominant halo mass quenching for the 
central galaxies. 


5 COMPARING THE HYBRID AND HALO 
QUENCHING MODELS 

In this section we perform a robust comparison between the 
two quenching models in two ways, an internal one based 
on Bayesian Information Criterion (BIC) described in Sec¬ 
tion 5.1, and an external one based on cross-validation (Sec¬ 
tion 5.2), which is motivated by the quenching impact on 
the average halo mass of the massive blue galaxies quanti¬ 
tatively explained in Section 4.3. 


5.1 Internal Model Gomparison: Bayesian 
Information Criterion 


In Bayesian applications, pairwise comparisons between 
models M[ and M 2 are often based on the Bayes factor 
6 i 2 , which is dehned as the ratio of the posterior odds, 
P(Mi\x)/P(M 2 \x), to the prior odds, n{Mi)ln{M 2 ). In our case, 
the bayes factor is 

„ _ B(^(,aio|x) ^(l2hybrid) 

“ Blghybndlx) ^ ’ 

SO that a 612 above unity indicates the data favor halo 
quenching and a 612 below points to hybrid quenching. How¬ 
ever, in most practical settings (as is the case here) the prior 
odds are hard to set precisely, and model selection based on 
BIC is widely employed as a rough equivalent to selection 
based on Bayes factors. The BIC (a.k.a., “Schwarz informa¬ 
tion criterion”), is defined as 

BIC =-21ni;ma,+ /tlnu, (34) 


where In Xmax is the maximum likelihood value, k is the num¬ 
ber of parameters, and n is the number of data points. Kass 
& Raftery (1995) argued that in the limit of large n (n=439 
in our analyses). 


-2 In 612 - (BIChalo - BlChybrid) 
-2 In 6|2 


(35) 


i.e., ABIC = BIChalo - BlChybrid can be viewed as a rough ap¬ 
proximation to - 21 nBi 2 , so that ABIC<0(< -10) indicates 
that 2haio is favored (strongly) and ABIC>0(> 10) points 
(strongly) to Qbybnd- 

The ABIC between the two quenching models is -35.8, 
which corresponds to an asymptotic value of lnBi2=17.9 ac¬ 
cording to equation 35. Therefore, based on the BIC test, 
the clustering and the g-g lensing measurements of the red 
and blue galaxies strongly favor the halo quenching model 
against the hybrid quenching model, and the halo mass is 
the more statistically dominant driver of galaxy quenching 
than stellar mass. 

The two quenching models are non-nested models with 
the same k and n, so the second term of equation 34 that 
penalizes model complexities is the same in both quenching 
models. The BIC test is then equivalent to the alternative 
Akaike information criterion that is based on relative likeli¬ 
hoods, or a simple Ay^ test (i.e., Ay^=35.8). These tests all 
point to the halo mass as the main driver of quenching. 


5.2 External Model Comparison: Halo Masses of 
Massive Blue Centrals 

The discussion in Section 4.3 points to a potentially 
smoking-gun test of the two quenching models, by compar¬ 
ing the host halo mass of the massive red and blue central 
galaxies predicted from the two best-fit quenching models, 
to other mass measurements for observed groups/clusters 
with red and blue centrals within the same redshift range. 
Unfortunately, clusters with blue centrals are systematically 
under-selected by most of the photometric cluster finders 
based on matching to the red sequence, while spectroscopic 
group catalogues constructed from friends-of-friends algo¬ 
rithms do not have large enough volume for finding many 
massive clusters. 

Recently, Mandelbaum (2015) constructed a sample 
of locally brightest galaxies (LBGs) from the SDSS main 
galaxy sample, by adopting a set of isolation criteria care¬ 
fully calibrated against semi-analytic mock galaxy cata¬ 
logues to minimize the satellite contamination rate (Wang 
2015). The resulting LBG sample is thus a subset of all 
massive central galaxies, but with excellent purity of cen¬ 
tral galaxy membership and zero bias against blue colour. 
Therefore, the LBGs are ideal for our purpose of measur¬ 
ing the segregation in halo mass between the red and blue 
centrals. 

Mandelbaum (2015) measured the average host halo 
mass of the LBGs directly by fitting an NFW density pro- 
Hle (after projection to 2D) to the weak lensing signals mea¬ 
sured below 1/z^'Mpc. Fig. 13 compares the host halo mass 
measured as a function of LBG stellar mass (data points 
with errorbars) to that predicted by the best-fit halo (solid) 
and hybrid (dashed) quenching models. The errorbars on the 
LBG measurements are the l-cr uncertainties on the mean 
halo mass, derived from 1000 bootstrap-resampled datasets. 
The coloured bands about the solid curves are the uncer¬ 
tainties on the mean halo mass predicted from the 68 % con- 
hdence regions. To avoid clutter, we do not show the un¬ 
certainties on the hybrid quenching predictions, which are 
comparable to the halo quenching uncertainties. The average 
halo mass predicted by the halo quenching model is in excel¬ 
lent agreement with the measurements from the LBG sam¬ 
ple, while the hybrid quenching model, as expected, grossly 
over-predicts the halo mass for the massive blue galaxies. 

The difference in the host halo mass between red and 
blue centrals, as predicted by the hybrid quenching model, 
can be understood simply as the outcome of a larger differ¬ 
ential growth between dark matter and stellar mass in quies¬ 
cent systems than in star-forming ones. More specifically, in 
quenched systems the dark matter halos usually continued 
to grow after the shutdown of stellar mass growth, while in 
star-forming systems the two often grew in sync, creating a 
bimodality in host halo mass between two colours without 
invoking halo quenching. However, this differential growth 
effect causes at most a factor of two difference in the av¬ 
erage host halo masses, too small to explain the factor of 
several difference observed at the high mass end (Quadri 
2015). 

The LBG experiment in Fig. 13 further demonstrates 
that the halo quenching model, employing halo mass as the 
driver for galaxy quenching, is superior to the hybrid quench¬ 
ing model, which relies on stellar mass to quench central 
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Figure 13. The average weak lensing mass of the host halos of red and blue central galaxies as functions of stellar mass. Solid and 
dashed thick curves are the predictions from the halo and hybrid quenching models, respectively. The colored bands show the the l-cr 
uncertainties about the two mean halo-to-stellar mass relations from the fiducial halo quenching model. The circles and squares with 
errorbars are the directly measured halo masses for the red and blue locally brightest galaxies (LBGs), respectively. The halo quenching 
model provides excellent agreement with both the red and blue LBG measurements, while the hybrid quenching model predicts a much 
higher halo mass for the blue LBGs. 


galaxies. As we explained in Section 4.3, the deficiency of 
the hybrid model in describing the signals of the massive 
blue galaxies is intrinsic to the stellar mass quenching mech¬ 
anism, which fails to explain the rare occurrence of blue cen¬ 
trals in massive clusters. The combined evidence from the 
BIG model comparison and the LBG experiment strongly 
suggests that the halo mass is the main driver for quenching 
the galaxies observed in SDSS. 


6 PHYSICAL IMPLICATIONS OF THE 
CONSTRAINTS ON HALO QUENCHING 
MODEL 

With the halo quenching model being established as the 
more viable scenario, we now focus back on the physical 
implications of our constraints on halo quenching. 

6.1 Uniform Characteristic Halo Masses for 
Quenching Centrals and Satellites 

Although the halo quenching formula for centrals and satel¬ 
lites are decoupled in the analysis, our fiducial constraint 
nonetheless recovers two very similar characteristic halo 


masses and for both species at around 1.5 x 

It is very tempting to associate this uniform 
quenching mass scale for both central and satellites to Msi,ock) 
the critical halo mass responsible for the turning-on of shock 
heating. Analytical calculations and hydrodynamic simula¬ 
tions both favor a Mj^ock of ~ few times (Birn- 

boim & Dekel 2003; Keres et al. 2005; Dekel & Birnboim 
2006), providing one of the most plausible explanation for 
the similar values of and derived statistically in our 
analyses. 

Conservatively speaking, even if the similarity between 
our inferred characteristic halo masses and Ms^ock were coin¬ 
cidental, the consistency between and still indicates 
that the quenching of centrals and satellites are somewhat 
coupled, most likely driven by processes that are both tied 
to the potential well of the halos. For instance, the super- 
massive black holes (SMBHs) could provide the thermal or 
mechanical feedback required to stop the halo gas from cool¬ 
ing and feeding the satellites (Di Matteo et al. 2005; Cro¬ 
ton et al. 2006; Somerville et al. 2008), while regulating the 
growth of the central galaxies (Ferrarese & Merritt 2000; 
Gebhardt et al. 2000; Tremaine et al. 2002). Hopkins et al. 
(2007) suggested that the SMBH mass is largely determined 
by the depth of the potential well in the central regions of 
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Figure 14. Comparison between the color-segregated stellar mass functions predicted by the best-fit halo quenching model (thick solid 
curves with shaded l-cr uncertainty bands) and that observed in the SDSS (shaded histograms) at three different redshifts (from left to 
right: z=0.07, 0.13, and 0.19). Within each color, the total SMF is further decomposed into the central (dashed) and satellite (dotted) 
contributions. The excellent agreement between the predicted and observed SMFs for both the red and blue populations is not a result 
of fitting, but the consistency between galaxy abundance and their clustering and lensing signals predicted by our iHQD analysis using 
the halo quenching model within ACDM. Note that above 10'*/i“^Mo, the observed SMF in the lowest redshift bin (left; z=0.07) is subject 
to slight incompleteness while the other two are relatively complete, hence a slight deviation from our best-fit prediction. 


the system, which precedes the assembly of halo mass, i.e., 
the maximum circular velocity is already half the present- 
day value by the time the halo has accreted only two per 
cent of its final mass (Bosch et al. 2014). 


6.2 Implications for Satellite Quenching and 
Galactic Conformity 

The halo quenching of satellites has a slower transition 
across than that of the central galaxies (left panel 

of Fig. 12). This rnles out the possibility that halo quench¬ 
ing does not distinguish between centrals and satellites. As 
mentioned in the introduction, even in the hot halo quench¬ 
ing scenario where gas cooling of centrals and satellites were 
equally inhibited, the satellites might experience significant 
delays in their quenching, due to a shorter exposure to the 
hot halo and/or a spell of star formation from pockets of cold 
gas they carried across the virial radius of the larger, hotter 
halo (Simha et al. 2009; Wetzel et al. 2012, 2013; van den 
Bosch et al. 2008). Additionally, recent observations suggest 
that other processes like pre-processing during infall (Haines 
et al. 2015), strangulation (Peng et al. 2015), and ram pres¬ 
sure stripping (Muzzin et al. 2014) are all at play, contribut¬ 
ing to the satellite quenching trend with halo mass. The in¬ 
efficiency of satellite quenching is also seen in dwarf galaxies 
below the stellar mass scale we probed here (Wheeler et al. 
2014). 

Another interesting aspect of the halo quenching sce¬ 
nario is that it may help explain galactic conformity, i.e., the 
observed correlation between colours of the central galaxies 
and their surrounding satellites (Weinmann et al. 2006; Kno- 


bel et al. 2015), because the quiescent pairs of centrals and 
satellites are quenched by the common halos they share. 
However, this halo quenching-induced conformity only oc¬ 
curs among central-satellite pairs at fixed M, of the cen¬ 
trals. To explain the galactic conformity observed at fixed 
Ml,, there either has to be a substantial scatter between ob¬ 
served Ml, and true Mi,, or a secondary process that couples 
the quenching of centrals and satellites within the same halo. 
For instance, halos formed earlier with higher concentration 
may be more likely to host quenched pairs of centrals and 
satellites than their younger and less concentrated coun¬ 
terparts at the same Mi, (Paranjape et al. 2015). Galactic 
conformity does not appear when the centrals and satellites 
were quenched independently, e.g., in the hybrid quenching 
scenario. 

The combination of this intra-halo conformity and the 
correlation between clustering bias and halo mass, could po¬ 
tentially explain the inter-halo conformity observed in the 
galaxy marked correlation statistics (Skibba et al. 2006; 
Cohn & White 2014) and hinted by galaxy pairs in the local 
volume (z<0.03; see Kauffmann et al. 2013), although a sec¬ 
ondary quenching induced by either formation time or halo 
concentration at fixed Mi, may be required (Paranjape et al. 
2015). We will explore the conformity prediction of the halo 
quenching model in the upcoming third paper of this series. 


6.3 Colour-segregated Stellar Mass Functions of 
Central and Satellite Galaxies 

As another consistency-check of our analysis. Fig. 14 shows 
the underlying SMFs of the red and blue galaxies pre- 
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Figure 15. Red (left) and blue (right) satellite stellar mass functions conditioned at six different host halo masses, predicted by the 
best-fit halo quenching model. For each conditional SMF curve, the shaded band indicates the l-cr uncertainty and the accompanying 
dashed curve indicate the conditional SMF of the combined red and blue satellites. The two inset panels show similar sets of curves as 
the main panel, but normalized by the mass of the corresponding host halos. 


dieted by our best-fit halo quenching model at three red- 
shifts (z=0.07, 0.13, 0.19). In each panel, the shaded bands 
are the l-cr uncertainties on the predicted SMFs, and the 
shaded histograms show the observed SMFs of each colour, 
i.e., direct galaxy number counts at each redshift without 
the 1/Vmax weighting. The SMF of each colour is further de¬ 
composed into contributions from the central (dashed) and 
satellite (dotted) galaxies. The observed SMF at z=0.07 (left 
panel) is more incomplete at the high M, end due to photo¬ 
metric confusions about bright sources in SDSS (see Paper 
I for details), therefore lying further below our predictions 
than that at the two higher redshifts. Since the observed 
SMFs are not used as input data to the constraints, the 
excellent agreement between our predictions and the direct 
number counts on the high stellar mass end is very encour¬ 
aging — it demonstrates the great consistency between the 
three key observables (i.e., galaxy clustering, g-g lensing, and 
the SMFs) measured for the two coloured populations and 
that predicted by the best-fit halo quenching model within 
the iHOD framework. Furthermore, the successful recovery 
of the red and blue SMFs means that the halo quenching 
model naturally recovers the stellar mass quenching trend 
observed in PIO. 

The best-fit halo quenching also provides concrete pre¬ 
diction for the conditional SMF of the red and blue satellite 
galaxies, shown on the left and right panels in Fig. 15, re¬ 
spectively. The conditional SMF is defined as the average 
number of satellites per dex in stellar mass at fixed halo 
mass, (dA^sat/dlgM,|M*), the integration of which over M, 
gives the commonly-used satellite HOD, i.e., the average 
number of galaxies per halo above some stellar mass limit 
(A'sat(fH,>Ml‘"')|M;,>. In each panel, the solid curves are the 


red/blue satellite SMFs within halos of six different masses 
( with their l-cr uncertainties indicated by the shaded 

bands. The dashed curves are the same in both panels, show¬ 
ing the two-colour combined satellite SMFs. Each inset panel 
shows the same set of curves as in the main panel, but 
each normalized by the corresponding halo mass. Overall, 
the satellite population above is dominated by the red 
galaxies, and the number of blue satellite galaxies per halo 
mass decreases with increasing halo mass due to progres¬ 
sively stronger halo quenching effect, while the total number 
of satellites per halo mass remains roughly constant. 


7 COMPARISON TO PREVIOUS WORKS 

Our quenching model is fundamentally different from pre¬ 
vious studies of the link between galaxy colours and the 
underlying dark matter halos. Here we compare our best-fit 
quenching model to the two main alternative methods. One 
is the separate red and blue galaxy modelling using tradi¬ 
tional HOD methods (Section 7.1), and the other is based 
on a modified abundance matching scheme, i.e., the age¬ 
matching model (Section 7.2). We summarize the compari¬ 
son between our result and the previous studies in Fig. 16, 
which zooms in on the stellar mass range of Fig. 13 that 
has the maximum model discriminating power. In addition 
to the LEG weak lensing masses shown in Fig. 13, we also 
include the average halo mass of the red and blue centrals 
measured from satellite kinematics by More et al. (2011). 
We note that although the various constraints and measure¬ 
ments shown in Fig. 13 assumed slightly different cosmolo¬ 
gies, the uncertainties due to cosmology are usually much 
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Figure 16. Comparison of various predicted and measured {Mi,\M,) relations for the red and blue galaxies. Thick solid curves with shaded 
bands are the predictions from our best-fit halo quenching model, while thick short-dashed curves are predicted from the age-matching 
model in Hearin et al. (2014). Long-dashed curves are the traditional HOD predictions from the Rodn'guez-Puebla et al. (2015) method 
which derived the red and blue relations from SDSS separately, while dot-dashed curves are from Tinker et al. (2013) for COSMOS 
galaxies at z>0.3. Filled and open symbols with errorbars are the average halo masses measured from the satellite kinematics of galaxy 
groups (More et al. 2011) and the weak lensing of LBGs (Mandelbaum 2015), respectively. 


smaller than the statistical errors, and the strong bimodal¬ 
ity (or the lack thereof) in the host halo mass between red 
and blue is independent of any changes in cosmology. We 
will come back to Fig. 16 frequently and discuss individual 
comparisons in detail below. 

7.1 Comparison to Traditional HOD Models 

The most straightforward way to model the red and blue 
split of galaxy observables traditionally is to infer the HODs 
of the overall and the red galaxies first, and subtract the two 
to derive the HOD of the blue. This approach guarantees the 
consistency between the three sets of HODs, but lacks the 
flexibility in the treatment of the red fraction for describing 
the full ranges of behaviours seen in the data (Zehavi et al. 
2005; Zehavi et al. 2011). A more comprehensive method is 


to treat the two colours separately, by prescribing indepen¬ 
dent HODs for the two and an ID overall quenched fraction 
as a function of halo mass, as done recently in Tinker et al. 
(2013) and Rodriguez-Puebla et al. (2015). 

As mentioned in the introduction, there are two main 
differences between our approach and the methods of Tin¬ 
ker et al. (2013) and Rodriguez-Puebla et al. (2015). Firstly, 
our quenching analysis employs only four more parameters 
to explain the split into red and blue galaxies, while the tra¬ 
ditional methods require doubling of the number of param¬ 
eters used for the overall population (e.g., 23 parameters in 
the Rodriguez-Puebla et al. 2015 analysis, and 27 in Tinker 
et al. 2013). Our four parameters are also more physically 
meaningful because they can be directly related to the aver¬ 
age quenching action, as we discussed in Section 6. Secondly, 
our quenching model describes the bimodality of galaxy oc- 
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cupation statistics in a mathematically consistent manner — 
the overall galaxy HOD is recovered by summing the inferred 
red and blue HODs. The other two methods do not benefit 
from this consistency, making the connection between the 
colour-segregated HODs and the overall HOD hard to in¬ 
terpret. Naively one might think that our model is a subset 
of the traditional HOD models that appear more flexible by 
fitting red and blue separately. However, since the combi¬ 
nation of the two separate coloured SHMRs derived in the 
traditional methods usually does not fit the overall galaxy 
population while our iHOD quenching models do, the two 
methods are fundamentally different descriptions of the red 
and blue galaxy populations. 

It is worth nothing that the constraints drawn from tra¬ 
ditional HOD modelling are unaffected by its incapability 
of including stellar mass quenching, which is subdominant 
compared to halo mass quenching. However, as we discussed 
in Section 6.2, the halo quenching of satellites has a much 
slower transition across the critical halo mass than that of 
centrals. Therefore, the lack of separate treatments for the 
red fractions in centrals and satellites remains an important 
issue in those traditional HOD models. 

Using the SMFs, galaxy clustering, and g-g lensing 
within the COSMOS survey. Tinker et al. (2013) derived 
the SHMRs of active and quiescent galaxies over the redshift 
range between 0.2 and 1.0. The active/quiescent classifica¬ 
tion is based on their separation on the optical vs. near-IR 
colour space, which is better at distinguishing dusty and SF 
objects than using the optical colours alone. Employing the 
global HOD framework of Leauthaud et al. (2012) they ap¬ 
plied independent central galaxy SHMR and satellite HODs 
to the two colours, and assumed a non-parametric form for 
the red fraction as a function of halo mass using as a spline- 
interpolated function through five pivotal halo masses. For 
the lowest redshift bin in their analysis (z~0.36), the derived 
red and blue SHMRs are very similar below M*~1O*^/!“*M0, 
but strongly diverge on the high Mi, end, with the red cen¬ 
trals having a lower average M, than the blue ones at fixed 
Ml,. This divergence is equivalent to having a stellar mass 
quenching at fixed halo mass, albeit in the opposite direc¬ 
tion of the stellar mass quenching trend observed in PIO. 
The {Mi,\M,) relation reveals a similar trend that our halo 
quenching model predicts, with the red central galaxies re¬ 
siding in more massive halos than the blue centrals, but the 
predicted difference between the red and blue amplitudes is 
much smaller (dot-dashed curves). 

Rodrfguez-Puebla et al. (2015) derived the separate 
SHMRS of red and blue galaxies using the combination 
of galaxy clustering and SMFs in SDSS. Instead of using 
the g-g lensing as an input, they employed the SMFs mea¬ 
sured for the centrals and satellites separately within each 
colour, using the SDSS group catalogue constructed by Yang 
et al. (2012). Unlike Tinker et al. (2013), they assumed a 
parametric form for the red fraction as a function of halo 
mass. The long-dashed curves in Fig. 16 are the re- 


^ The parameterisation of the IHDD framework in Paper I is also 
heavily based on the Leauthaud et al. (2012) framework, but is 
fundamentally different in the treatment of sample completeness 
and signal calculation; see Paper I for a detailed comparison be¬ 
tween the two frameworks. 


lations inferred by Rodrfguez-Puebla et al. (2015), showing 
good agreement with our predictions as well as the measure¬ 
ments from satellite kinematics and LEG weak lensing. The 
higher amplitude of {Mi,\M,) on M,>2 x 1O*'Y^M0 is driven 
by the extrapolation of the parametric relation from lower 
M, rather than data. Although Rodrfguez-Puebla et al.’s and 
our analyses make use of the same galaxies in SDSS, the two 
constraints are derived using independent methods and dif¬ 
ferent measurements. Hence, the good degree of consistency 
between these results is non-trivial. 

7.2 Comparison to the Age-Matching Model 

An alternative way to describe the colour dependence of 
galaxy clustering and g-g lensing is to extend the SHAM 
model to allow a secondary matching between galaxy 
colour and subhalo formation time, i.e., the so-called “age¬ 
matching” model (Hearin & Watson 2013). In particular, 
after the usual mapping between subhalo mass (using Vpeak 
as a proxy; see Reddick et al. 2013) and galaxy stellar mass, 
the age-matching method rank-orders the characteristic red¬ 
shift Zstarve of tliose subhalos at fixed M, and matches the 
galaxies hosted by older subhalos to redder colours, while 
maintaining the observed colour distribution of galaxies at 
that M,. For most of the centrals (below 10**/i“^Mo), Zstarve 
is equivalent to the formation redshift of the subhalos Zfornu 
but at very high M, it is dominated by Zchan the first epoch 
at which halo mass exceeds IO'^/i^'Mq. The age-matching 
method roughly reproduces the colour and stellar mass de¬ 
pendences of the clustering and g-g lensing signals (Hearin 
et al. 2014). 

In the language of statistical quenching, the age¬ 
matching method describes the red fraction as 

f,t{M„ Ml,) = X (36) 

where ni(zstarve|44*) is determined by the matching between 
^starve and colour at fixed M,, and ZstarvcCfH/i) is the formation 
time vs. halo mass relation. Therefore, the age-matching pro¬ 
cess to first order assumes a stellar mass quenching, as the 
colour-matching step is done in bins of M, regardless of halo 
mass (the first term on the RHS of equation 36), while the 
secondary quenching is via formation time (the second term 
on the RHS of equation 36). The combined quenching effect 
is best illustrated in Fig. 17, where we show the distribution 
of red fraction on the Mf-Mi, diagram for centrals (left) and 
satellites (right), calculated from the age-matching mock 
catalogue generated by Hearin et al. (2014) (with the origi¬ 
nal NYU-VAGC stellar mass and halo virial mass). For the 
centrals, the stellar mass quenching along the vertical axis 
dominates, albeit in a discrete fashion due to the binning 
artefact. At fixed stellar mass, since halo age is a decreas¬ 
ing function of halo mass (i.e., ^Mi,<Q-, c.f., figure 12 

in Wechsler et al. 2002), the secondary quenching direction 
of the age-matching method is a reversed halo quenching, 
with bluer centrals occupying younger, thus more massive 
halos. Therefore, according to the discussion in Section 4.3, 
we anticipate the “age-quenching” to predict rela¬ 

tions that are similar to the hybrid quenching model, i.e., a 
weak segregation in the host halo mass between the red and 
blue centrals. 

The magenta and cyan dotted curves in Fig. 16 indicate 
the {Mi,\M,) relations of the red and blue galaxies, measured 
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Figure 17. Similar to Fig. 4, but for the quenching within the age-matching model. The colour map shows the red galaxy fraction as a 
function of both stellar mass and halo mass for the central (left) and satellite (right) galaxies in the age-matching mock catalogue. Unlike 
other figures in this paper, we use the VAGC stellar mass and the halo virial mass used in the original age-matching catalogue. Clearly 
the centrals exhit strong stellar mass quenching trend, with a secondary formation-time (i.e., age) quenching trend at fixed M,, while 
the satellites are dominated by stellar mass qunching. The reversed halo quenching due to age-matching explains the overprediction of 
host halo mass for blue centrals. 


from the age-matching mock catalogue produced by Hearin 
et al. (2014) (with both M, and M* converted to our mass 
units and definitions). Compared to the halo mass measured 
from satellite kinematics and LBG weak lensing, the age¬ 
matching mock heavily over-predicts the amplitude of the 
relation for the blue galaxies, making the red and blue cen¬ 
trals occupy halos of similar mass, despite the systematic 
difference in their formation times. The tiny difference in the 
average halo mass between the two colours (below 0.2 dex 
over all mass scales) also switches sign across M,~10*'/r“^Mo, 
where the indicator for Zstarve switches from Zfom, to Zchar- In 
particular, below 10**/i“^Mo, age-matching predicts that the 
blue centrals live in more massive halos than the red cen¬ 
trals, due to the anti-correlation between Zfonn and halo mass 
— at any given M, the younger halos that are assigned bluer 
centrals are also more massive. However, this increase of halo 
mass with bluer colour of the centrals is in the opposite di¬ 
rection compared to the observations. Above 10'*/!“^Mo, the 
characteristic redshift Zstarve is dominated by Zchar, which is 
more positively correlated with halo mass and assigns redder 
centrals to more massive systems. The overall disagreement 
between the {Mi,\M,) predicted by the age-matching method 
and that measured from satellite kinematics and LBG weak 
lensing indicates that the stellar mass quenching assumed in 
age-matching is not adequate, and the secondary formation¬ 
time quenching at fixed stellar mass is strongly disfavored 
by the observations. 

Another difference between the age-matching model 
and the quenching models considered in this paper is that, 
by choosing formation time as a quenching indicator, the 
age-matching model exhibits the maximum level of galactic 
assembly bias (Zentner et al. 2014), which is absent in our 


quenching models, by construction. However, while the “halo 
assembly bias”, namely, the dependence of halo properties 
on the formation history, is clearly detected in cosmological 
simulations (Sheth & Tormen 2004; Gao et al. 2005; Wech- 
sler et al. 2006; Harker et al. 2006; Zhu et al. 2006; Jing 
et al. 2007; Hahn et al. 2007; Groft et al. 2012), whether it 
left a significant imprint on the observed galaxy properties is 
still in debate (Berlind et al. 2006; Blanton & Berlind 2007; 
Wang et al. 2013; Miyatake et al. 2015; Lin et al. 2015). 
In mock galaxy catalogues constructed from semi-analytical 
models and hydrodynamic simulations, galaxy clustering ex¬ 
hibits an assembly bias of at most ~ 10 per cent (Yoo et al. 
2006; Groton et al. 2007; Zu et al. 2008; Mehta 2014).^ The 
great success of the halo quenching model in quantitatively 
explaining the clustering and weak lensing of the red and 
blue galaxies, while simultaneously recovering their respec¬ 
tive SMFs and average halo masses, strongly suggests that 
the halo quenching is the dominant process in shaping the 
distribution of galaxy colours observed in SDSS, and that 
any impact of the galactic assembly bias should be a sec¬ 
ondary effect, in the form of, e.g., a formation-time (or con¬ 
centration) quenching at fixed halo mass proposed by Paran- 
jape et al. (2015). 

Finally, the physical interpretation of the halo quench¬ 
ing model (as discussed in Section 6.1) relates the quenching 
of galaxies to the capability of the host halos to either heat 


^ Some hydrodynamic simulations predict much higher assembly 
bias effect on clustering (~20%), which cannot be captured by 
any current abundance matching models (Chaves-Montero et al. 
2015). 
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the incoming gas or keep the hot gas from cooling, and the 
sharp transition of this capability across some critical halo 
mass is the key to explain the strong bimodality in galaxy 
colours. This physical picture is fundamentally different from 
the physical motivation of age-matching, in which the galaxy 
quenching is strictly tied to the dark matter accretion his¬ 
tory of halos at fixed M,. However, the average halo accre¬ 
tion history is a smooth function of cosmic time, therefore 
showing no bimodality in the formation time of halos (Zhao 
et al. 2009). In addition, the connection between dark matter 
accretion and SFR is very complex. Using a suite of high- 
resolution hydrodynamic simulations, Faucher-Giguere et al. 
(2011) showed that the cold gas accretion rate, which is more 
directly related to star formation, is in general not a simple 
universal factor of the dark matter accretion rate, and that 
baryonic feedback can cause SFRs to deviate significantly 
from the external gas accretion rates. 


8 CONCLUSIONS 

We develop a novel method to identify the dominant driver 
of galaxy quenching in the local Universe, using the galaxy 
clustering and g-g lensing of red and blue galaxies observed 
in SDSS. The method extends the powerful iHOD frame¬ 
work developed in Zu & Mandelbaum (2015) by introducing 
two quenching models: 1) a halo quenching model in which 
the average probability of a galaxy being quenched depends 
solely on the main halo mass, but in separate manners for 
centrals and satellites; and 2) a hybrid quenching model in 
which the quenching probability of all galaxies depends on 
their stellar mass, with the satellite having an extra depen¬ 
dence on the host halo mass. 

The two quenching models predict distinctive 2D distri¬ 
butions of red galaxy fractions on the M,-M* plane, result¬ 
ing in different patterns through which red and blue galaxies 
populate dark matter halos. We then predict the clustering 
and g-g lensing signals of the red and blue galaxies from 
these two galaxy occupation patterns using the iHOD frame¬ 
work and compare them to the measurements from SDSS. 
Most importantly, the flexibility of iHOD allows us to include 
~80% more galaxies in the analysis than the traditional HOD 
method, greatly enhancing our capability of statistically dis¬ 
tinguishing the two quenching models. 

We find that the halo quenching model provides better 
descriptions of the bimodality in the clustering and lensing of 
observed galaxies than the hybrid quenching models, mainly 
due to the significantly improved fit to massive blue galax¬ 
ies. We further identify that the average host halo mass of 
the massive blue centrals provides the most discriminating 
power in testing viable quenching models — models with 
halo mass quenching generally predict a much stronger seg¬ 
regation in the average host halo mass ((M;,|M,)) between 
the red and blue than the ones without (Fig. 13). 

Therefore, by comparing the predicted by var¬ 

ious quenching models, including the age-matching model, 
to that measured directly from the satellite kinematics of 
galaxy clusters and the weak lensing of locally brightest 
galaxies, we conhrm that the best-fit halo quenching model 
provides excellent agreement with the two observational 
measurements, while models that rely on stellar mass (e.g., 
the hybrid quenching model and the age-matching method) 


fail to predict the halo mass of the blue central galax¬ 
ies (Fig. 16). Furthermore, the formation time-quenching at 
hxed M, prescribed in the age-matching method creates a 
reversed halo quenching trend, therefore placing blue and 
red centrals of the same M, into higher and lower mass ha¬ 
los, respectively. This trend is strongly disfavored by the 
observations where at any given M, redder centrals on aver¬ 
age occupy more massive halos. Our findings indicate that 
any viable abundance matching scheme for assigning galaxy 
colours has to reproduce the observed strong bimodality in 
host halo mass between red and blue centrals, especially at 
the high mass end. 

The derived characteristic halo masses of the central 
and satellite quenching have very similar values around 
1.5 X 10'^/!“*Mo, suggesting that a uniform halo quenching 
process is operating on both the centrals and satellites. The 
derived characteristic halo mass can be interpreted by the 
canonical halo quenching theory, which predicts a critical 
halo mass of Msi,ock~10*^/!“'MQ. Above this critical mass, the 
virial shock is able to prevent star formation by heating the 
infalling gas to high temperatures, whereas below Mshock the 
halo quenching rapidly turns off, creating strong bimodal¬ 
ity in both the colour and the spatial distribution of galax¬ 
ies. The pace at which the halo mass quenching operates, 
however, appears to be faster for the centrals than for the 
satellite galaxies, which are quenched in a more delayed and 
prolonged fashion. 

In the future, we anticipate that the iHOD halo quench¬ 
ing model, which accurately explains the spatial clustering, 
g-g lensing, and SMFs of the red and blue galaxies, will 
provide an important baseline model for explaining an even 
wider range of observed galaxy properties with ever-growing 
precision. In the near term, i.e., the upcoming Paper III in 
this series, we plan to generate realistic colour-segregated 
galaxy mock catalogues using the constraints inferred from 
the halo quenching model analysis in this paper, and make 
a comprehensive comparison with the observed galaxies to 
look for potential signatures of any secondary quenching pro¬ 
cesses, e.g., due to formation time (Paranjape et al. 2015) 
and galaxy compactness (Woo et al. 2015) at fixed halo mass. 
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